Exploration of the potential common pathogenic mechanisms in COVID-19 and silicosis by using bioinformatics and system biology.

Silicosis is an occupational lung disease that is common worldwide. In recent years, coronavirus disease 2019 (COVID-19) has provided daunting challenges to public healthcare systems globally. Although multiple studies have shown a close link between COVID-19 and other respiratory diseases, the inter-relational mechanisms between COVID-19 and silicosis remain unclear. This study aimed to explore the shared molecular mechanisms and drug targets of COVID-19 and silicosis. Gene expression profiling identified four modules that were most closely associated with both diseases. Furthermore, we performed functional analysis and constructed a protein-protein interaction network. Seven hub genes (budding uninhibited by benzimidazoles 1 [BUB1], protein regulator of cytokinesis 1 [PRC1], kinesin family member C1 [KIFC1], ribonucleotide reductase regulatory subunit M2 [RRM2], cyclin-dependent kinase inhibitor 3 [CDKN3], Cyclin B2 [CCNB2], and minichromosome maintenance complex component 6 [MCM6]) were involved in the interaction between COVID-19 and silicosis. We investigated how diverse microRNAs and transcription factors regulate these seven genes. Subsequently, the correlation between the hub genes and infiltrating immune cells was explored. Further in-depth analyses were performed based on single-cell transcriptomic data from COVID-19, and the expression of hub-shared genes was characterized and located in multiple cell clusters. Finally, molecular docking results reveal small molecular compounds that may improve COVID-19 and silicosis. The current study reveals the common pathogenesis of COVID-19 and silicosis, which may provide a novel reference for further research.

Bioinformatics and system biology approach to identify the influences among COVID-19, ARDS and sepsis.

Background Severe coronavirus disease 2019 (COVID -19) has led to severe pneumonia or acute respiratory distress syndrome (ARDS) worldwide. we have noted that many critically ill patients with COVID-19 present with typical sepsis-related clinical manifestations, including multiple organ dysfunction syndrome, coagulopathy, and septic shock. The molecular mechanisms that underlie COVID-19, ARDS and sepsis are not well understood. The objectives of this study were to analyze potential molecular mechanisms and identify potential drugs for the treatment of COVID-19, ARDS and sepsis using bioinformatics and a systems biology approach. Methods Three RNA-seq datasets (GSE171110, GSE76293 and GSE137342) from Gene Expression Omnibus (GEO) were employed to detect mutual differentially expressed genes (DEGs) for the patients with the COVID-19, ARDS and sepsis for functional enrichment, pathway analysis, and candidate drugs analysis. Results We obtained 110 common DEGs among COVID-19, ARDS and sepsis. ARG1, FCGR1A, MPO, and TLR5 are the most influential hub genes. The infection and immune-related pathways and functions are the main pathways and molecular functions of these three diseases. FOXC1, YY1, GATA2, FOXL, STAT1 and STAT3 are important TFs for COVID-19. mir-335-5p, miR-335-5p and hsa-mir-26a-5p were associated with COVID-19. Finally, the hub genes retrieved from the DSigDB database indicate multiple drug molecules and drug-targets interaction. Conclusion We performed a functional analysis under ontology terms and pathway analysis and found some common associations among COVID-19, ARDS and sepsis. Transcription factors-genes interaction, protein-drug interactions, and DEGs-miRNAs coregulatory network with common DEGs were also identified on the datasets. We believe that the candidate drugs obtained in this study may contribute to the effective treatment of COVID-19.

Using bioinformatics and systems biology methods to identify the mechanism of interaction between COVID-19 and nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) is considered a risk factor for severe COVID-19, but the mechanism remains unknown. This study used bioinformatics to help define the relationship between these diseases. The GSE147507 (COVID-19), GSE126848 (NAFLD), and GSE63067 (NAFLD-2) datasets were screened using the Gene Expression Omnibus. Common differentially expressed genes were then identified using a Venn diagram. Gene ontology analysis and KEGG pathway enrichment were performed on the differentially expressed genes. A protein-protein interaction network was also constructed using the STRING platform, and key genes were identified using the Cytoscape plugin. GES63067 was selected for validation of the results. Analysis of ferroptosis gene expression during the development of the 2 diseases and prediction of their upstream miRNAs and lncRNAs. In addition, transcription factors (TFs) and miRNAs related to key genes were identified. Effective drugs that act on target genes were found in the DSigDB. The GSE147507 and GSE126848 datasets were crossed to obtain 28 co-regulated genes, 22 gene ontology terms, 3 KEGG pathways, and 10 key genes. NAFLD may affect COVID-19 progression through immune function and inflammatory signaling pathways. CYBB was predicted to be a differential ferroptosis gene associated with 2 diseases, and the CYBB-hsa-miR-196a/b-5p-TUG1 regulatory axis was identified. TF-gene interactions and TF-miRNA coregulatory network were constructed successfully. A total of 10 drugs, (such as Eckol, sulfinpyrazone, and phenylbutazone) were considered as target drugs for Patients with COVID-19 and NAFLD. This study identified key gene and defined molecular mechanisms associated with the progression of COVID-19 and NAFLD. COVID-19 and NAFLD progression may regulate ferroptosis through the CYBB-hsa-miR-196a/b-5p-TUG1 axis. This study provides additional drug options for the treatment of COVID-19 combined with NAFLD disease.

IMC-MDA: Prediction of miRNA-disease association based on induction matrix completion.

To comprehend the etiology and pathogenesis of many illnesses, it is essential to identify disease-associated microRNAs (miRNAs). However, there are a number of challenges with current computational approaches, such as the lack of "negative samples", that is, confirmed irrelevant miRNA-disease pairs, and the poor performance in terms of predicting miRNAs related with "isolated diseases", i.e. illnesses with no known associated miRNAs, which presents the need for novel computational methods. In this study, for the purpose of predicting the connection between disease and miRNA, an inductive matrix completion model was designed, referred to as IMC-MDA. In the model of IMC-MDA, for each miRNA-disease pair, the predicted marks are calculated by combining the known miRNA-disease connection with the integrated disease similarities and miRNA similarities. Based on LOOCV, IMC-MDA had an AUC of 0.8034, which shows better performance than previous methods. Furthermore, experiments have validated the prediction of disease-related miRNAs for three major human diseases: colon cancer, kidney cancer, and lung cancer.

Multi-pathogen based chimeric vaccine to fight against COVID-19 and concomitant coinfections.

BACKGROUND: COVID-19 has proved to be a fatal disease of the year 2020, due to which thousands of people globally have lost their lives, and still, the infection cases are at a high rate. Experimental studies suggested that SARS-CoV-2 interacts with various microorganisms, and this coinfection is accountable for the augmentation of infection severity. METHODS AND RESULTS: In this study, we have designed a multi-pathogen vaccine by involving the immunogenic proteins from S. pneumonia, H. influenza, and M. tuberculosis, as they are dominantly associated with SARS-CoV-2. A total of 8 antigenic protein sequences were selected to predict B-cell, HTL, and CTL epitopes restricted to the most prevalent HLA alleles. The selected epitopes were antigenic, non-allergenic, and non-toxic and were linked with adjuvant and linkers to make the vaccine protein more immunogenic, stable, and flexible. The tertiary structure, Ramachandran plot, and discontinuous B-cell epitopes were predicted. Docking and MD simulation study has shown efficient binding of the chimeric vaccine with the TLR4 receptor. CONCLUSION: The in silico immune simulation analysis has shown a high level of cytokines and IgG after a three-dose injection. Hence, this strategy could be a better way to decrease the disease's severity and could be used as a weapon to prevent this pandemic.

Identifying key genes related to inflammasome in severe COVID-19 patients based on a joint model with random forest and artificial neural network.

Background: The coronavirus disease 2019 (COVID-19) has been spreading astonishingly and caused catastrophic losses worldwide. The high mortality of severe COVID-19 patients is an serious problem that needs to be solved urgently. However, the biomarkers and fundamental pathological mechanisms of severe COVID-19 are poorly understood. The aims of this study was to explore key genes related to inflammasome in severe COVID-19 and their potential molecular mechanisms using random forest and artificial neural network modeling. Methods: Differentially expressed genes (DEGs) in severe COVID-19 were screened from GSE151764 and GSE183533 via comprehensive transcriptome Meta-analysis. Protein-protein interaction (PPI) networks and functional analyses were conducted to identify molecular mechanisms related to DEGs or DEGs associated with inflammasome (IADEGs), respectively. Five the most important IADEGs in severe COVID-19 were explored using random forest. Then, we put these five IADEGs into an artificial neural network to construct a novel diagnostic model for severe COVID-19 and verified its diagnostic efficacy in GSE205099. Results: Using combining P value < 0.05, we obtained 192 DEGs, 40 of which are IADEGs. The GO enrichment analysis results indicated that 192 DEGs were mainly involved in T cell activation, MHC protein complex and immune receptor activity. The KEGG enrichment analysis results indicated that 192 GEGs were mainly involved in Th17 cell differentiation, IL-17 signaling pathway, mTOR signaling pathway and NOD-like receptor signaling pathway. In addition, the top GO terms of 40 IADEGs were involved in T cell activation, immune response-activating signal transduction, external side of plasma membrane and phosphatase binding. The KEGG enrichment analysis results indicated that IADEGs were mainly involved in FoxO signaling pathway, Toll-like receptor, JAK-STAT signaling pathway and Apoptosis. Then, five important IADEGs (AXL, MKI67, CDKN3, BCL2 and PTGS2) for severe COVID-19 were screened by random forest analysis. By building an artificial neural network model, we found that the AUC values of 5 important IADEGs were 0.972 and 0.844 in the train group (GSE151764 and GSE183533) and test group (GSE205099), respectively. Conclusion: The five genes related to inflammasome, including AXL, MKI67, CDKN3, BCL2 and PTGS2, are important for severe COVID-19 patients, and these molecules are related to the activation of NLRP3 inflammasome. Furthermore, AXL, MKI67, CDKN3, BCL2 and PTGS2 as a marker combination could be used as potential markers to identify severe COVID-19 patients.

Identification of biomarkers, pathways, and potential therapeutic targets for heart failure using next-generation sequencing data and bioinformatics analysis.

BACKGROUND: Heart failure (HF) is the most common cardiovascular diseases and the leading cause of cardiovascular diseases related deaths. Increasing molecular targets have been discovered for HF prognosis and therapy. However, there is still an urgent need to identify novel biomarkers. Therefore, we evaluated biomarkers that might aid the diagnosis and treatment of HF. METHODS: We searched next-generation sequencing (NGS) dataset (GSE161472) and identified differentially expressed genes (DEGs) by comparing 47 HF samples and 37 normal control samples using limma in R package. Gene ontology (GO) and pathway enrichment analyses of the DEGs were performed using the g: Profiler database. The protein-protein interaction (PPI) network was plotted with Human Integrated Protein-Protein Interaction rEference (HiPPIE) and visualized using Cytoscape. Module analysis of the PPI network was done using PEWCC1. Then, miRNA-hub gene regulatory network and TF-hub gene regulatory network were constructed by Cytoscape software. Finally, we performed receiver operating characteristic (ROC) curve analysis to predict the diagnostic effectiveness of the hub genes. RESULTS: A total of 930 DEGs, 464 upregulated genes and 466 downregulated genes, were identified in HF. GO and REACTOME pathway enrichment results showed that DEGs mainly enriched in localization, small molecule metabolic process, SARS-CoV infections, and the citric acid tricarboxylic acid (TCA) cycle and respiratory electron transport. After combining the results of the PPI network miRNA-hub gene regulatory network and TF-hub gene regulatory network, 10 hub genes were selected, including heat shock protein 90 alpha family class A member 1 (HSP90AA1), arrestin beta 2 (ARRB2), myosin heavy chain 9 (MYH9), heat shock protein 90 alpha family class B member 1 (HSP90AB1), filamin A (FLNA), epidermal growth factor receptor (EGFR), phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1), cullin 4A (CUL4A), YEATS domain containing 4 (YEATS4), and lysine acetyltransferase 2B (KAT2B). CONCLUSIONS: This discovery-driven study might be useful to provide a novel insight into the diagnosis and treatment of HF. However, more experiments are needed in the future to investigate the functional roles of these genes in HF.

Analysis of Persian Bioinformatics Research with Topic Modeling.

Purpose: As a scientific field, bioinformatics has drawn remarkable attention from various fields, such as information technology, mathematics, and modern biological sciences, in recent years. The topic models originating from the field of natural language processing have become the focus of attention with the rapid accumulation of biological datasets. Thus, this research is aimed at modeling the topic content of the bioinformatics literature presented by Iranian researchers in the Scopus Citation Database. Methodology. This research was a descriptive-exploratory study, and the studied population included 3899 papers indexed in the Scopus database, which had been indexed in this database until March 9, 2022. The topic modeling was then performed on the abstracts and titles of the papers. A combination of LDA and TF-IDF was utilized for topic modeling. Findings. The data analysis with topic modeling resulted in identifying seven main topics "Molecular Modeling," "Gene Expression," "Biomarker," "Coronavirus," "Immunoinformatics," "Cancer Bioinformatics," and "Systems Biology." Moreover, "Systems Biology" and "Coronavirus" had the largest and smallest clusters, respectively. Conclusion: The present investigation demonstrated an acceptable performance for the LDA algorithm in classifying the topics included in this field. The extracted topic clusters indicated excellent consistency and topic connection with each other.

Harnessing Deep Learning for Omics in an Era of COVID-19.

Omics data are multidimensional, heterogeneous, and high throughput. Robust computational methods and machine learning (ML)-based models offer new prospects to accelerate the data-to-knowledge trajectory. Deep learning (DL) is a powerful subset of ML inspired by brain structure and has created unprecedented momentum in bioinformatics and computational biology research. This article provides an overview of the current DL models applied to multi-omics data for both the beginner and the expert user. Additionally, COVID-19 will continue to impact planetary health as a pandemic and an endemic disease, with genomic and multi-omic pathophysiology. DL offers, therefore, new ways of harnessing systems biology research on COVID-19 diagnostics and therapeutics. Herein, we discuss, first, the statistical ML algorithms and essential deep architectures. Then, we review DL applications in multi-omics data analysis and their intersection with COVID-19. Finally, challenges and several promising directions are highlighted going forward in the current era of COVID-19.

MSClustering: A Cytoscape Tool for Multi-Level Clustering of Biological Networks.

MSClustering is an efficient software package for visualizing and analyzing complex networks in Cytoscape. Based on the distance matrix of a network that it takes as input, MSClustering automatically displays the minimum span clustering (MSC) of the network at various characteristic levels. To produce a view of the overall network structure, the app then organizes the multi-level results into an MSC tree. Here, we demonstrate the package's phylogenetic applications in studying the evolutionary relationships of complex systems, including 63 beta coronaviruses and 197 GPCRs. The validity of MSClustering for large systems has been verified by its clustering of 3481 enzymes. Through an experimental comparison, we show that MSClustering outperforms five different state-of-the-art methods in the efficiency and reliability of their clustering.

Galaxy Training: A powerful framework for teaching!

There is an ongoing explosion of scientific datasets being generated, brought on by recent technological advances in many areas of the natural sciences. As a result, the life sciences have become increasingly computational in nature, and bioinformatics has taken on a central role in research studies. However, basic computational skills, data analysis, and stewardship are still rarely taught in life science educational programs, resulting in a skills gap in many of the researchers tasked with analysing these big datasets. In order to address this skills gap and empower researchers to perform their own data analyses, the Galaxy Training Network (GTN) has previously developed the Galaxy Training Platform (https://training.galaxyproject.org), an open access, community-driven framework for the collection of FAIR (Findable, Accessible, Interoperable, Reusable) training materials for data analysis utilizing the user-friendly Galaxy framework as its primary data analysis platform. Since its inception, this training platform has thrived, with the number of tutorials and contributors growing rapidly, and the range of topics extending beyond life sciences to include topics such as climatology, cheminformatics, and machine learning. While initially aimed at supporting researchers directly, the GTN framework has proven to be an invaluable resource for educators as well. We have focused our efforts in recent years on adding increased support for this growing community of instructors. New features have been added to facilitate the use of the materials in a classroom setting, simplifying the contribution flow for new materials, and have added a set of train-the-trainer lessons. Here, we present the latest developments in the GTN project, aimed at facilitating the use of the Galaxy Training materials by educators, and its usage in different learning environments.

Identification of key differentially expressed genes in SARS-CoV-2 using RNA-seq analysis with a systems biology approach.

COVID-19 is associated with dysregulation of several genes and signaling pathways. Based on the importance of expression profiling in identification of the pathogenesis of COVID-19 and proposing novel therapies for this disorder, we have employed an in silico approach to find differentially expressed genes between COVID-19 patients and healthy controls and their relevance with cellular functions and signaling pathways. We obtained 630 DEmRNAs, including 486 down-regulated DEGs (such as CCL3 and RSAD2) and 144 up-regulated DEGs (such as RHO and IQCA1L), and 15 DElncRNAs, including 9 down-regulated DElncRNAs (such as PELATON and LINC01506) and 6 up-regulated DElncRNAs (such as AJUBA-DT and FALEC). The PPI network of DEGs showed the presence of a number immune-related genes such as those coding for HLA molecules and interferon regulatory factors. Taken together, these results highlight the importance of immune-related genes and pathways in the pathogenesis of COVID-19 and suggest novel targets for treatment of this disorder.

Exploration of the Shared Molecular Mechanisms between COVID-19 and Neurodegenerative Diseases through Bioinformatic Analysis.

The COVID-19 pandemic has caused millions of deaths and remains a major public health burden worldwide. Previous studies found that a large number of COVID-19 patients and survivors developed neurological symptoms and might be at high risk of neurodegenerative diseases, such as Alzheimer's disease (AD) and Parkinson's disease (PD). We aimed to explore the shared pathways between COVID-19, AD, and PD by using bioinformatic analysis to reveal potential mechanisms, which may explain the neurological symptoms and degeneration of brain that occur in COVID-19 patients, and to provide early intervention. In this study, gene expression datasets of the frontal cortex were employed to detect common differentially expressed genes (DEGs) of COVID-19, AD, and PD. A total of 52 common DEGs were then examined using functional annotation, protein-protein interaction (PPI) construction, candidate drug identification, and regulatory network analysis. We found that the involvement of the synaptic vesicle cycle and down-regulation of synapses were shared by these three diseases, suggesting that synaptic dysfunction might contribute to the onset and progress of neurodegenerative diseases caused by COVID-19. Five hub genes and one key module were obtained from the PPI network. Moreover, 5 drugs and 42 transcription factors (TFs) were also identified on the datasets. In conclusion, the results of our study provide new insights and directions for follow-up studies of the relationship between COVID-19 and neurodegenerative diseases. The hub genes and potential drugs we identified may provide promising treatment strategies to prevent COVID-19 patients from developing these disorders.

The molecular mechanism of cardiac injury in SARS-CoV-2 infection: Focus on mitochondrial dysfunction.

BACKGROUND: Coronavirus disease 2019(COVID-19) caused a large number of infections worldwide. Although some patients recovered from the disease, some of the other problems that accompanied it, such as cardiac injury, could affect the patient's subsequent quality of life and prognosis. OBJECTIVES: To clarify the molecular mechanism of cardiac injury in SARS-CoV-2 Infection. METHODS: The RNA-Seq dataset (GSE184715) comparing expression profiling of Mock human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) and SARS-CoV-2-infected hiPSC-CMs was downloaded from Gene Expression Omnibus (GEO). Differentially expressed genes(DEGs) were performed by the R software. Degs were analyzed by enrichment analysis to clarify the affected pathways. Hub genes were screened out by a PPI network constructed from Degs. Finally, Connectivity Map was used to screen for the treatment of COVID-19 induced cardiac injury. RESULTS: 2705 differentially expressed genes were identified. Enrichment analysis confirmed that mitochondrial dysfunction was caused by SARS-CoV-2, meanwhile, cardiac muscle contraction was suppressed and NF-κB was activated. Based on the PPI network, 15 hub genes were identified. These 15 down-regulated hub genes were mainly involved in the reduced activity of complexes in the mitochondrial respiratory chain associated with mitochondrial dysfunction. Moreover, 5 candidate drugs were identified to treat cardiac injury. CONCLUSION: In conclusion, SARS-CoV-2 infection of cardiomyocytes causes mitochondrial dysfunction, including reduced mitochondrial respiratory chain complex activity and decreased ATP synthesis, leading to cardiomyocyte apoptosis, while the activated NF-κB also induced cytokine storms, ultimately resulting in cardiac injury.

Immunoinformatics-aided design of a new multi-epitope vaccine adjuvanted with domain 4 of pneumolysin against Streptococcus pneumoniae strains.

BACKGROUND: Streptococcus pneumoniae (Pneumococcus) has remained a leading cause of fatal infections such as pneumonia, meningitis, and sepsis. Moreover, this pathogen plays a major role in bacterial co-infection in patients with life-threatening respiratory virus diseases such as influenza and COVID-19. High morbidity and mortality in over one million cases, especially in very young children and the elderly, are the main motivations for pneumococcal vaccine development. Due to the limitations of the currently marketed polysaccharide-based vaccines, non-serotype-specific protein-based vaccines have received wide research interest in recent years. One step further is to identify high antigenic regions within multiple highly-conserved proteins in order to develop peptide vaccines that can affect various stages of pneumococcal infection, providing broader serotype coverage and more effective protection. In this study, immunoinformatics tools were used to design an effective multi-epitope vaccine in order to elicit neutralizing antibodies against multiple strains of pneumococcus. RESULTS: The B- and T-cell epitopes from highly protective antigens PspA (clades 1-5) and PhtD were predicted and immunodominant peptides were linked to each other with proper linkers. The domain 4 of Ply, as a potential TLR4 agonist adjuvant candidate, was attached to the end of the construct to enhance the immunogenicity of the epitope vaccine. The evaluation of the physicochemical and immunological properties showed that the final construct was stable, soluble, antigenic, and non-allergenic. Furthermore, the protein was found to be acidic and hydrophilic in nature. The protein 3D-structure was built and refined, and the Ramachandran plot, ProSA-web, ERRAT, and Verify3D validated the quality of the final model. Molecular docking analysis showed that the designed construct via Ply domain 4 had a strong interaction with TLR4. The structural stability of the docked complex was confirmed by molecular dynamics. Finally, codon optimization was performed for gene expression in E. coli, followed by in silico cloning in the pET28a(+) vector. CONCLUSION: The computational analysis of the construct showed acceptable results, however, the suggested vaccine needs to be experimentally verified in laboratory to ensure its safety and immunogenicity.

Identification of a key glioblastoma candidate gene, FUBP3, based on weighted gene co-expression network analysis.

BACKGROUND: Glioblastoma multiforme (GBM) is the most common aggressive malignant brain tumor. However, the molecular mechanism of glioblastoma formation is still poorly understood. To identify candidate genes that may be connected to glioma growth and development, weighted gene co-expression network analysis (WGCNA) was performed to construct a gene co-expression network between gene sets and clinical characteristics. We also explored the function of the key candidate gene. METHODS: Two GBM datasets were selected from GEO Datasets. The R language was used to identify differentially expressed genes. WGCNA was performed to construct a gene co-expression network in the GEO glioblastoma samples. A custom Venn diagram website was used to find the intersecting genes. The GEPIA website was applied for survival analysis to determine the significant gene, FUBP3. OS, DSS, and PFI analyses, based on the UCSC Cancer Genomics Browser, were performed to verify the significance of FUBP3. Immunohistochemistry was performed to evaluate the expression of FUBP3 in glioblastoma and adjacent normal tissue. KEGG and GO enrichment analyses were used to reveal possible functions of FUBP3. Microenvironment analysis was used to explore the relationship between FUBP3 and immune infiltration. Immunohistochemistry was performed to verify the results of the microenvironment analysis. RESULTS: GSE70231 and GSE108474 were selected from GEO Datasets, then 715 and 694 differentially expressed genes (DEGs) from GSE70231 and GSE108474, respectively, were identified. We then performed weighted gene co-expression network analysis (WGCNA) and identified the most downregulated gene modules of GSE70231 and GSE108474, and 659 and 3915 module genes from GSE70231 and GSE108474, respectively, were selected. Five intersection genes (FUBP3, DAD1, CLIC1, ABR, and DNM1) were calculated by Venn diagram. FUBP3 was then identified as the only significant gene by survival analysis using the GEPIA website. OS, DSS, and PFI analyses verified the significance of FUBP3. Immunohistochemical analysis revealed FUBP3 expression in GBM and adjacent normal tissue. KEGG and GO analyses uncovered the possible function of FUBP3 in GBM. Tumor microenvironment analysis showed that FUBP3 may be connected to immune infiltration, and immunohistochemistry identified a positive correlation between immune cells (CD4 + T cells, CD8 + T cells, and macrophages) and FUBP3. CONCLUSION: FUBP3 is associated with immune surveillance in GBM, indicating that it has a great impact on GBM development and progression. Therefore, interventions involving FUBP3 and its regulatory pathway may be a new approach for GBM treatment.

Identification of differentially expressed genes and pathways in kidney of ANCA-associated vasculitis by integrated bioinformatics analysis.

The antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV) is a systematic of relatively rare autoimmune diseases with unknown cause. Kidney involvement is one of the most common clinical manifestations, and the degree of renal damage is closely associated with the development and prognosis of AAV. In this study, we utilized the Robust Rank Aggreg (RRA) method in R to integrate GSE104948, GSE104954, GSE108109, GSE108112, and GSE108113 profile datasets loaded from Gene Expression Omnibus (GEO) database and identified a set of differentially expressed genes (DEGs) in kidney between AAV patients and living donors. Then, the results of gene ontology (GO) functional annotation showed that immunity and metabolism involved process of AAV both in glomerulus and tubulointerstitial. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that following pathways, such as complement and coagulation cascades pathway; Staphylococcus aureus infection; disease-COVID-19; and systemic lupus erythematosus (SLE) pathway play a crucial role in AAV. Next, the results analyzed by protein-protein interaction (PPI) network and Cytoscape software exhibited the hub genes ALB, TYROBP, and CYBB existed in both glomerular and tubulointerstitial compartments datasets. Finally, KEGG analysis using genes of two most important modules also further validated complement and coagulation cascades pathway and S. aureus infection existed both in glomerulus and tubulointerstitial compartments datasets. In conclusion, this study identified key genes and pathways involved in kidney of AAV, which was benefit to further uncover the mechanisms underlying the development and progress of AAV, biomarkers, and potential therapeutic targets as well.

Design of a multi-epitope-based vaccine consisted of immunodominant epitopes of structural proteins of SARS-CoV-2 using immunoinformatics approach.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has shown rapid global spread and has resulted in a significant death toll worldwide. In this study, we aimed to design a multi-epitope vaccine against SARS-CoV-2 based on structural proteins S, M, N, and E. We identified B- and T-cell epitopes and then the antigenicity, toxicity, allergenicity, and similarity of predicted epitopes were analyzed. T-cell epitopes were docked with corresponding HLA alleles. Consequently, the selected T- and B-cell epitopes were included in the final construct. All selected epitopes were connected with different linkers and flagellin and pan-HLA DR binding epitopes (PADRE) as an adjuvant were used in the vaccine construct. Furthermore, molecular docking was used to evaluate the complex between the final vaccine construct and two alleles, HLA-A*02:01 and HLA-DRB1*01:01. Finally, codons were optimized for in silico cloning into pET28a(+) vector using SnapGene. The final vaccine construct comprised 11 CTL, HTL, and B-cell epitopes corresponding to 394 amino acid residues. In silico evaluation showed that the designed vaccine might potentially promote an immune response. Further in vivo preclinical and clinical testing is required to determine the safety and efficacy of the designed vaccine.

Bioinformatic analysis and preliminary validation of potential therapeutic targets for COVID-19 infection in asthma patients.

BACKGROUND: Severe acute respiratory syndrome coronavirus 2 causes coronavirus disease 19 (COVID-19). The number of confirmed cases of COVID-19 is also rapidly increasing worldwide, posing a significant challenge to human safety. Asthma is a risk factor for COVID-19, but the underlying molecular mechanisms of the asthma-COVID-19 interaction remain unclear. METHODS: We used transcriptome analysis to discover molecular biomarkers common to asthma and COVID-19. Gene Expression Omnibus database RNA-seq datasets (GSE195599 and GSE196822) were used to identify differentially expressed genes (DEGs) in asthma and COVID-19 patients. After intersecting the differentially expressed mRNAs, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to identify the common pathogenic molecular mechanism. Bioinformatic methods were used to construct protein-protein interaction (PPI) networks and identify key genes from the networks. An online database was used to predict interactions between transcription factors and key genes. The differentially expressed long noncoding RNAs (lncRNAs) in the GSE195599 and GSE196822 datasets were intersected to construct a competing endogenous RNA (ceRNA) regulatory network. Interaction networks were constructed for key genes with RNA-binding proteins (RBPs) and oxidative stress-related proteins. The diagnostic efficacy of key genes in COVID-19 was verified with the GSE171110 dataset. The differential expression of key genes in asthma was verified with the GSE69683 dataset. An asthma cell model was established with interleukins (IL-4, IL-13 and IL-17A) and transfected with siRNA-CXCR1. The role of CXCR1 in asthma development was preliminarily confirmed. RESULTS: By intersecting the differentially expressed genes for COVID-19 and asthma, 393 common DEGs were obtained. GO and KEGG enrichment analyses of the DEGs showed that they mainly affected inflammation-, cytokine- and immune-related functions and inflammation-related signaling pathways. By analyzing the PPI network, we obtained 10 key genes: TLR4, TLR2, MMP9, EGF, HCK, FCGR2A, SELP, NFKBIA, CXCR1, and SELL. By intersecting the differentially expressed lncRNAs for COVID-19 and asthma, 13 common differentially expressed lncRNAs were obtained. LncRNAs that regulated microRNAs (miRNAs) were mainly concentrated in intercellular signal transduction, apoptosis, immunity and other related functional pathways. The ceRNA network suggested that there were a variety of regulatory miRNAs and lncRNAs upstream of the key genes. The key genes could also bind a variety of RBPs and oxidative stress-related genes. The key genes also had good diagnostic value in the verification set. In the validation set, the expression of key genes was statistically significant in both the COVID-19 group and the asthma group compared with the healthy control group. CXCR1 expression was upregulated in asthma cell models, and interference with CXCR1 expression significantly reduced cell viability. CONCLUSIONS: Key genes may become diagnostic and predictive biomarkers of outcomes in COVID-19 and asthma. Video Abstract.